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Abstract 



We investigate the effects of randomness in a strongly correlated electron 
model in one-dimension at half-filling. The ground state correlation functions 
are exactly written by products of 3x3 transfer matrices and are evaluated 
numerically. The correlation lengths depend on randomness when the inter- 
action is effectively weak. On the contrary, they are completely insensitive to 
randomness when the interaction is effectively strong. 
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The behavior of electrons in the presence of randomness has attracted a lot of attentions 
as one of the most fundamental problems in condensed matter physics Q. In the absence 
of interaction, the scaling theory gives us a criterion whether states are localized or not 
||. It was shown rigorously that all the states are localized for wide class of models in one 
dimension [[3] . In two dimensions, it is believed that all the states are localized. Randomness 
induces a metal-insulator transition in three dimensions. So far, a lot of theoretical and 
experimental works have been reported. However, the validity to describe the experiments 
by non-interacting models is an open question, since Coulomb interaction between electrons 
is always present. 

In one dimension, some interacting models without randomness can be solved exactly by 
the Bethe ansatz technique or bosonization 0] and the properties have been investigated. 
However, approaches of including randomness to such models seem to be hopeless. In 
the presence of randomness and without interaction, exact results were obtained on the 
localization of eigenstates ||. In this way, the individual investigations of models with 
interaction or randomness have been succeeded considerablely. However, it is extremely hard 
task to take into account them simultaneously and we have no reliable method to analize. 
Although a few results by the perturbation method [0J or by bosonization |J are known, 
even the qualitative understandings are far from satisfactory. Numerical investigations have 
difficulty due to the restriction of the system size. Furthermore, one needs enormous amount 
of CPU time for averaging over samples to obtain enough accuracy. 

In this paper, we study a special model at half-filling to avoid the numerical difficulties 
mentioned above to investigate the effect of randomness in the strongly correlated electron 
model. The lattice structure is shown in Fig. [I] and the Hamiltonian with the open boundary 
condition is given by 

L 

H = V E { E [( " PiaPia - tipUia ~ Up\+lJia + h.C.) + Vfnfa ]+<+ < +l J\V , (1) 
t=1,I i=l 

where a unit cell is labeled by i. Here pi a is the annihilation operator with spin a =f, J. at 
site i with no interaction (a p-site) which have at most two electrons with opposite spins, 
di a is the the annihilation operator at a site with infinitely large on-site Coulomb repulsion 
(a d-site) which can have at most one electron, and nf a (a = p, d) is the electron number 
operator. The projection operator which represents the infinitely large on-site Coulomb 
repulsion on d-sites is V = 0^(1 — n it n ii)- We denote the on-site potentials for rf-sites by 
V s. For simplicity we parametrize t's and V d 's by positive A's as 

U = Xi (2) 
yf = + 2. (3) 

Then the on-site potentials at p-sites are set to be zero except at the boundaries. We shall 
take A's to be independent random variable. The advantages of the model (0) is that the 
exact and unique ground state (at half-filling) is explicitly written as 

*o)=^n n {pi + pIi« + vt) |o). (4) 

»=i<t=T4 
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Without randomness, namely when A's are uniform, the exact ground state in a restricted 
parameter space was obtained || by following the construction introduced by Brandt and 
Giesekus M. The correlation functions and the momentum distribution functions were 
calculated exactly The correlation functions are exactly represented by products of 

the 3x3 transfer matrices P,|T0f. Therefore, the correlation functions and the correlation 
lengths can be numerically obtained for considerably long chains even for random chains. 

When there is neither interaction nor randomness, the ground state is a band insulator. 
With interaction and without randomness, the ground state is also an insulator but totally 
different type due to the existence of a spin gap |TTJ . Then the model enables us to investigate 
the effect of randomness to an insulating state (where the insulating behavior is due to the 
strong correlation). The extension of the method to other parametrizations of the model 
and other models of similar type is straightforward. 

We take uniform randomness with width W 



A-y<A,(=^)<A + _ 



(5) 



The probability density function for W < 2 A is 



p(x = V* 



W^8(2-x) 

otherwise 



for -2 (A + f ) 2 + 2 < x < -2 (A - f ) 2 + 2 



(6) 



and for W > 2A 



p{x = V?) 



HVfey for-2(A + f) 2 + 2<x<-2(A-f) 2 + 2 

, 2 



2== for -2 (A - f) + 2 < x < 2 

l(2-x) V 2 / — 



W y/2(2-x) 

otherwise 



(7) 



For all W, average is 



V? = -2 (A 2 + ^)+2 



The difference of the on-site potentials between p- and d-sites depends both on A and W. 

Without randomness, the spin, density, singlet-pair, and (cl^Cj^ correlation functions 
decay exponentially f|]||. This suggests the existence of a finite excitation gap above the 
ground state and it was numerically confirmed ||11|| . Of course, on d-sites Coulomb interaction 
is always infinity. However, by choosing on-site potential V d , the model interpolates between 
the following two limits: (i) A — > oo. In this limit V d — > — oo and (n d ) — > 1. Since each 
d-site is occupied by one electron, they forbid to have additional electrons. In this sense the 
effective interaction is strong, (ii) A — ► 0. In this limit V d — > 2. Since the hopping matrix 
elements between p- and rf-sites, A, is infinitesimal comparing with V d , one has (n d ) — > 0. 
Thus the effective interaction is weak since no electron is in ci-sites. 

The occupation and the correlation functions are exactly written M 
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where 



"1" 




"1" 
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, F = 





1 








(11) 



Here (9 a 's are the number, spin, creation (anihilation) of singlet-pair, or creation (ambula- 
tion) operators, and a, (3 = p or d. The matrices T's, M's, and O's are the corresponding 
transfer matrices is given by 



T 



M h = — 



or 
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(12) 



A^ + A,- Xf 



A? 
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for the spin correlation functions 
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for the density correlation functions 



(15) 
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Ma = 1, 
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for the singlet-pair correlation functions 



(16) 



The correlation lengths of the correlation functions between p-sites and between d-sites are 
the same up to order 0(1/L), since only the matrices at i and j sites are different in the 
representation (|10D. Due to the same reason, the correlation lengths of the spin and the 
singlet-pair correlation functions are the same up to order 0(1/ L). For a fixed set {Aj}, we 
numerically evaluate the quantities 

i L B +N 

— T (nd 

i=L B 



1 



L B +N 



m=j—i 



i=L B 



(17) 



where N is the number of sites which are used for the averaging in a sample and Lb is the 
number of sites which is ignored to exclude contributions form the boundary. We choose 
L = 10000, Lb = 2500, and N = 5000. The occupations of d-sites are shown in Fig. [| 
The sizes of the error bars are smaller than those of the plotted points. Note that (nf) 
= 2 — (nf), since the system is half-filled. We confirmed that the correlation functions decay 

exponentially. The correlation lengths are given from (OfOj) m _^_. oc exp [— m/t^o], where 



O = S for the spin and O = c for the correlation function (cj^c ^). The correlation lengths 



The estimates of £ c and 



are estimated by least squre fit for the values log 

£s are shown in Figs. |] and [|, respectively. The size of the error bars is smaller than that 
of the plotted points. 

The behavior of the correlation lengths depends on the occupation of d-sites, namely 
the effective interaction. For the parameter regime A < 1 where the effective interaction 
is weak, the correlation lengths become short as W increases. This behavior seems to be 
similar to the non-interacting cases. For A ^> 1 where the effective interaction is strong, the 
correlation lengths are independent of the strength of the randomness. We obtained similar 
behaviors for the density and the singlet-pair correlation functions. 

For the non-interacting cases, spin degree of freedom has nothing to do with the prop- 
erties of the systems. For the interacting cases, on the other hand, spin degree of freedom 
plays an important role and the effects of randomness are likely to be different from those 
for the non-interacting cases. The ground state (|J) is given by superpositions of spin singlet 
states. Within the analysis of this model, the results suggest that the states where the 
effective interaction is strong have a tendency to have local nature and the overlappings do 
not contribute much to the expectation values of the correlations. Thus their properties are 
stable against randomness. 
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FIGURES 

FIG. 1. The lattice structure. An open circle denotes a p-site (with no interaction) and a solid 
circle denotes a d-site (with infinitely large on-site Coulomb repulsion). A line represents hopping 
of electrons. 



FIG. 2. The estimates of {nf) as functions of IF for A = 0.1, 0.2, 0.5, 1, and 10. 



FIG. 3. The estimates of £ c as functions of IF for A = 0.1, 0.2, 0.5, 1, and 10. 



FIG. 4. The estimates of £ 5 as functions of IF for A = 0.1, 0.2, 0.5, 1, and 10. 



